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For an Ising spin glass on a hierarchical lattice, we show that the energy barrier to be overcome 
during the flip of a domain of size L scales as L d_1 for all dimensions d. We do this by investigating 
appropriate lower bounds to the barrier energy, which can be evaluated using an algorithm that 
remains fast for large system sizes and dimensions. The asymptotic limit d — > oo is evaluated 
analytically. 
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Energy barriers determine the dynamics of glassy sys- 
tems that have a complex energy landscape with many 
metastable states. Typically, the fluctuations between 
free energy minima in these systems (either between dif- 
ferent realizations or in the same random system) scale 
with observation size L as L e . It is generally assumed 
that the free energy barriers encountered in moving from 
one minimum to another scale with observation size as 
Lr . However, there exist so far few numerical studies 
of the value of this exponent in spin glasses, probably 
because of the difficulty of the problem (in contrast to 
the many studies of 9). For directed polymers in ran- 
dom systems, the identity tp = 9 was demonstrated some 
time ago For Ising spin glasses, Fisher and Huse 
derived within the droplet picture the double inequality 
9 < tp < d — 1 ,2] for dimension d. The lower limit is 
due to the fact that a domain wall has to be introduced 
into the system if all its spins are to be flipped. However, 
the minimum domain wall energy scales as L . The up- 
per limit is obtained by moving a straight domain wall 
through the system. Since such a domain wall breaks 
bonds, its energy cannot be larger than ~ L d_1 . 
Experiments on two- Q and three-dimensional spin 
glasses as well as a numerical studies in two dimensions 
|5j point towards a value of tp close to or identical to d— 1. 
On the other hand, an equality ip = 9 is sometimes tac- 
itly assumed, as for instance in a recent publication on 
spin glass dynamics on the hierarchical lattice |fj , where 
the probability for a spin flip on the length scale L is 
chosen to be a function of the effective coupling strength 
on this scale, which increases as L e . 

In this paper, we consider the Ising spin glass on an 
hierarchical lattice and show in fact that ip = d — 1 in all 
dimensions. The Hamiltonian of the system is given by 
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where (ij) indicates the sum over all nearest-neighbor 
pairs, and the spins assume the values ±1. We will 
mainly consider a Gaussian distribution of couplings with 



FIG. 1: Iterative construction of the hierarchical lattice. At 
each step, each bond is replaced with 2 d_1 pairs of bonds, with 
a new lattice site in between. Different grey shades indicate 
sites added at different iteration steps. 



zero mean and unit width. The hierarchical lattice gives 
for the exponent 9 in dimension d = 2 the value —0.27 
7] which compares well with the value —.29 from re- 
cent numerical studies on the square lattice [2j while in 
three dimensions the hierarchical lattice gives 9 w 0.25 
while on a simple cubic lattice its value is close to 
0.20 0. Thus the hierarchical lattice provides reason- 
ably good estimates for the value of the exponent 9 (at 
least for low-dimensional systems) and it is our hope that 
it is equally useful for determining the value of ip. 

The problem of finding energy barriers in g lassy sys- 
tems numerically is usually NP-complete [l0(. We will 
therefore not attempt to calculate the barrier exactly, 
but we will rather place bounds on it. Since the upper 
bound ip = d — 1 is already known due to the above-cited 
argument by Fisher and Huse, we will show in the follow- 
ing that there exists a lower bound to the barrier energy 
that increases with system size as L d ~ x . 

A hierarchical lattice is constructed by starting with 
one bond connecting two sites. This bond is replaced 
with a unit consisting of 2 d ~ 1 pairs of bonds, with a new 
site between each pair. Each of the 2 d bonds is again 
replaced with a unit of 2 d_1 pairs of bonds, etc., leading 
to a lattice with 2 Id bond after / iterations. In Fig.^the 
first three steps of this process are illustrated for d = 2. 
Evaluating a thermodynamic quantity on the hierarchical 
lattice with / levels is equivalent to evaluating it on a d- 
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dimensional hypercubic lattice of linear size 2 1 using the 
Migdal-Kadanoff-approximation. Recently, also the dy- 
namics of spin gl asses have been studied on hierarchical 
lattices [a [llj, [l2j , although there is no simple relation to 
the dynamics on hypercubic lattices. In approxi- 
mations based on renormalization ideas were made. 

We focus on the energy barrier that has to be overcome 
when moving from the ground state to the lowest-energy 
configuration with a domain wall. This domain-wall con- 
figuration is obtained by flipping one of the two level-0 
spins and by determining the new ground state with this 
fixed new configuration of the level-0 spins. We allow 
only single-spin flips when moving from the initial to the 
final state. At zero temperature (the situation we are 
considering here), the free-energy barrier is identical to 
the energy barrier. As indicated above, our goal is to 
show that there exists a lower bound that scales as L^ 1 . 
For this purpose, we consider the neighborhood of the 
right-hand level-0 spin, the "corner spin" (see Fig. |5J). 
We focus our attention on its L d ~ x level-/ nearest neigh- 
bors and its (L/2) d ~ 1 level- (7—1) next-nearest neighbors. 
At the moment where the corner spin is flipped, the next- 
nearest neighbors are in a configuration C nnn . We do not 
know the configuration of these spins which is associated 
with the true barrier, but we know that the optimum 
spin-flip sequence which passes through the true barrier 
state must have one of the possible configurations C nnn of 
the next-nearest neighbor spins at the moment where the 
corner spin is flipped. Therefore we will later minimize 
our lower bound with respect to C nnn . We start from the 
configuration of lowest energy that can be obtained with 
a given configuration C nnn of the next nearest neighbors 
of the right-hand corner spin and with the two corner 
spins fixed at their initial configuration. Clearly, the en- 
ergy of this state is at least as high as the energy of the 
initial configuration. We next calculate the minimum en- 
ergy barrier that has to be overcome when the right-hand 
corner spin is flipped with the configuration of the next- 
nearest neighbors remaining fixed. This minimum energy 
barrier is obtained by first flipping a suitable selection of 
the nearest-neighbor spins of the corner spin, before the 
corner spin itself is flipped. The barrier state is the one 
immediately before or immediately after the corner spin 
is flipped, and it is reached from our initial state by spin 
flips each of which increases the energy. Minimizing the 
energy barrier (i.e. the energy difference between the 
barrier state and our initial state) with respect to C nnn 
gives a lower bound to the true barrier. This is because 
the energy of our initial configuration is at least as high 
as that of the true initial configuration and since the en- 
ergy of our barrier state cannot be larger than that of 
the true barrier state. 

In the following, we determine this lower bound to the 
barrier. The right-hand corner spin is connected to each 
of its next-nearest neighbors via 2 d ~ 1 intermediate (level- 
I) spins, each of which initially assume the orientation 




FIG. 2: (a) The right-hand corner site drawn large, (b) its 
nearest-neighbors drawn large, (c) its next-nearest neighbors 
drawn large. 



that has the lower energy. Our task now consists in find- 
ing those intermediate spins that have to be flipped be- 
fore the corner spin is flipped, such that the barrier en- 
ergy is as low as possible. Intermediate spins, for which 
the absolute value of the coupling to its left neighbor is 
stronger than the absolute value of the coupling with the 
corner spin, must not be flipped. This is because for a 
given configuration of its two neighboring spins, the en- 
ergy is always lower when the intermediate spin has the 
orientation that satisfies the left-hand bond and possibly 
frustrates the right-hand bond. Since the values of the 
couplings are assigned at random, only about half of the 
intermediate spins are candidates for being flipped before 
the corner spin. 

Evaluating all possible combinations of these interme- 
diate spins that might be flipped before the corner spin 
costs a computer time that increases exponentially with 
the number of these spins. We will therefore later make 
an approximation that underestimates the above-defined 
lower bound to the barrier. In order to define the quanti- 
ties we need, let us first consider a subunit of three spins 
(see Fig. The corner spin on the right, one of its 
next-nearest neighbors to the left of it, and the interme- 
diate spin sitting between the two and connected to both 
of them. The left-hand spin is fixed. The corner spin 
is first in its initial configuration. The intermediate spin 
has the configuration that minimizes the energy. This ini- 
tial energy is our reference energy, and we set it to zero. 
Now let be the energy of the three-spin unit when 
the intermediate spin is flipped, and let e' 2 ' be the en- 
ergy of the three-spin unit when the right-hand spin (the 
corner spin) is flipped without first flipping the interme- 
diate spin, and let be the energy of the three-spin unit 
when the right-hand spin is flipped and the intermediate 
spin is adjusted such that it minimizes the energy. If the 
left-hand bond is stronger than the right-hand bond, we 
have e' 2 ' = otherwise we have e' 2 ' > e^K 

Now let n count all those three-spin units for which the 
intermediate spin is being flipped before the corner spin, 
and let m count all those three-spin units for which the 
intermediate spin is not flipped before the corner spin. 
The energy of the system just before the corner spin is 
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FIG. 3: A subunit of three spins and the energies associated 
with the different spin configurations mentioned in the text. 
The figure is drawn for the case that the right-hand bond is 
stronger than the left-hand bond (\J2\ > |J1|). 




FIG. 4: Ef >u b/2 d ~ 1 as function of dimension d. The dashed 
horizontal line indicates the analytically determined asymp- 
totic value 0.3305 which is approached in the limit d — ► oo. 



flipped is 
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and the energy of the system right after the corner spin 
is flipped is 
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The lower bound we are looking for is then 
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Since S underestimates Eb , it is also a lower bound to 
the barrier, and we focus on it in the following. We now 
perform the minimization with respect to the configura- 



"bubble" consisting of the corner spin, a next-nearest 
neighbor and 2 d ~ 1 intermediate spins, the contribution 
to the above sum is minimized if the next-nearest neigh- 
bor is in the configuration that leads to the higher initial 
bubble energy. We are interested in the average of the 
above sum over many different systems. Since the con- 
tributions of the bubbles to this average are additive, it 
is sufficient to take the average over one bubble, and our 
lower bond S is then for large L simply the number of 
bubbles, (L/2) d_1 , times this average. The lower bound 
to the barrier is therefore for large L given by 



S = (L/2) d - 1 x E, 
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i.e., it scales as L d ~ x . For a Gaussian distribution of unit 
width of the couplings, the bubble average has in three 
dimensions the value E bub = 0.364, i.e. it is positive. It 
increases with increasing dimension and approaches for 
large dimensions the asymptotic value E bub — 0.3305 x 
2 d " 1 (see Fig. EI). 

For the limit d — > oo we can prove analytically that 
Etub is positive. For each three-spin unit, we call the 
weaker coupling Jl, and the stronger coupling J2. In 
the limit d — > oo, one quarter of all three-spin units 
are initially not frustrated and have the stronger bond 
on the left-hand side. They make together a contribu- 
tion 2 d ~ 3 (|Jl|) to E bub . One quarter of all three-spin 
units are initially not frustrated and have the weaker 
bond on the left-hand side. They make together a con- 
tribution 2 d ~ 3 (| J2|) to E bub . One quarter of all three- 
spin units are initially frustrated in the right-hand bond 
and make together a contribution — 2 d ~ 3 (|Jl|) to E bub . 
The last quarter of all three-spin units are initially frus- 
trated in the left-hand bond and make together a con- 
tribution 2 d ~ 3 (| J2| - 2 1 Jl\) to E bub . Together this gives 
E bub = 2 d ~ 2 (|J2| - |J1|), which is positive since |J2| > 
|J1| 
Ull 



For a Gaussian bond distribution of unit width, 
= 2^%i = 0.467 and \ J2\ = 4= = 1.128, leading to 

E bub = 0.3305 x 2 d ~ 1 , in agreement with our asymptotic 
limit in Fig. 0] 

Our lower bound S obtained by calculation (J2J) is neg- 
ative for d = 2 and is therefore useless in this case. The 
reason is that the value of is more often negative in 
d = 2 than in higher dimensions (because we start from 
the higher initial bubble energy) , leading to a large differ- 
ence between E B ^ and E^ and making the inequality in 
the second line of © worse than in higher dimensions. In 
order to obtain a positive E bub also in d = 2. we replace 



in the second line of © the term \{E B 



E^>) with 



the more general expression aE B ^ + (1 — a)E B ' 
is valid for all < a < 1. We then obtain 
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tion of the next-nearest neighbor spins, C nnn . For each E B and makes E bub indeed positive. Already for a 



Choosing a > 0.5 places more weight on the larger energy 
(i) 
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FIG. 5: The figure shows a plot of E bub /2 d - 1 for a J = 
±2' 1-d '' /2 bond distribution after one decimation step as func- 
tion of dimension d. 

0.51 we obtain a positive Eb u f,/2 ~ 0.0039, and for a = 
0.74, it is as large as 0.12. In d = 1, the barrier energy 
is twice the absolute value of the largest bond, which 
increases as \/ln L for a Gaussian distribution of bonds. 
We therefore obtain (apart from logarithmic corrections 
in d = 1) the general result that the barrier scales on the 
hierarchical lattice for all d as L d_1 , implying ip = d— 1. 

In the following, we argue that our results can be gen- 
eralized to other bond distributions. For some distribu- 
tions, our expression for can become negative. This 
happens for instance for the ±J model, where the fi- 
nal state of the bubble can be reached from the initial 
state without first increasing the energy (Eh u b = —1-5 
for d = 3). In this case we modify our calculation in 
the following way: we first perform a sufficient amount 
of decimation steps on the hierarchical lattice, until the 
distribution of the renormalized bonds is such that our 
expression S (Eq. (J2J) for the lower bound becomes posi- 
tive. It must eventually become positive since the asymp- 
totic bond distribution (if rescaled to unit width) is for 
all distributions with finite moments the same and is very 
close to a Gaussian distribution. Then our argument can 
be repeated on the coarse-grained lattice, giving a new 
lower bound to the barrier energy for the original lattice. 
This is because the energy of the system with a given 
configuration of those spins that survived the renormal- 
ization procedure can never be smaller than that of the 
renormalized system. It turned out that performing one 
decimation step is sufficient for the ±J model even in 
d = 3 in order to obtain a positive value of E^. Fig. [S] 
shows Ehub/^ -1 as function of d for the distribution of 
bonds that is obtained after one decimation step for bond 
values J = ±2 < - 1 - d ^ 2 . (The values are chose such that 
the width of the distribution after one decimation step 
is 1.) The data are almost indistinguishable from those 
of Fig. 01 This is to be expected in higher dimensions, 



since the bond distribution after one decimation step is 
a binomial distribution, which approaches a Gaussian in 
high dimensions. We therefore obtain even for the ±J 
system a lower bound to the barrier that scales as L d ~ l . 
However, we can expect strong finite-size effects for small 
system sizes. 

To conclude, we have shown that the energy barrier 
that has to be overcome when introducing a domain wall 
into an Ising spin glass on a hierarchical lattice scales in 
all dimensions as L d_1 , for all bond distributions with 
finite moments. It remains to be seen if these results 
can be generalized to conventional lattices. However, we 
find it encouraging that the experimental data for ip seem 
strongly to favor a value close to d— 1 rather than the "ri- 
val" value 9. Since experimental results on spin glasses are 
not always probing droplets on the length scales at which 
droplet scaling ideas can be expected to apply without 
the use of corrections to scaling ||| , it is unrealistic to ex- 
pect perfect agreement for the value of the exponent ip. 
A further complication is that experimental and numeri- 
cal data are often studied at temperatures quite close to 
the transition temperature T c , when crossover to critical 
behavior will also complicate the analysis |l3( ■ 
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